KEK Preprint 2000-118 
NWU-HEP 2000-01 
DPNU-00-35 
TUAT-HEP 2000-5 
OCU-HEP 2000-01 

Measurement of the cross-section and forward-backward 
charge asymmetry for the b and c-quark in e + e 

annihilation with inclusive muons at yfs = 58 GeV 

TOPAZ Collaboration 

Y. Inoue 1 , A. Miyamoto 2 , E. Nakano 1 , T. Takahashi 1 , T. Tauchi 2 , Y. Tcramoto 1 , K. Abe 3 , T. 
Abe 3 [] I. Adachi 2 , K. Adachi 4 , M. Aoki 3 , M. Aoki 5 , R. Enomoto 2 ^ K. Emi 2 , H. Fujii 2 , K. Fujii 2 , 
T. Fujii 69 , J. Fujimoto 2 , N. Fujiwara 4 , H. Hirano 5 , B. Howell 7 , H. Hayashii 4 , N. Iida 2 , H. Ikeda 2 , 

S. Itami 3 , R. Itoh 2 , H. Iwasaki 2 , M. Iwasaki 4 ^, R. Kajikawa 3 , K. Kaneyuki 5 , S. Kato 2 , S. 
Kawabata 2 , H. Kichimi 2 , M. Kobayashi 2 , D. Koltick 7 , I. Levine 7 , H. Mamada 8 , K. Miyabayashi 4 , 
K. Muramatsu 4 , K. Nagai 9 Q, K. Nakabayashi 3 , M. Nakamura 1 , S. Noguchi 4 , O. Nitoh 8 , A. 
Ochi 5 , F. Ochiai 10 , N. Ohishi 3 , Y. Ohnishi 2 , Y. Ohshima 5 , H. Okuno 2 , T. Okusawa 11 , E. 
Shibata 7 , A. Sugiyama 3 , H. Sugiyama 4 , S. Suzuki 3 , K. Takahashi 2 , T. Tanimori 5 |, M. Tomoto 3 , 
T. Tsukamoto 2 , T. Tsumura 8 , S. Uno 2 , Y. Watanabe 5 , A. Yamamoto 2 , and M. Yamauchi 2 

1 Institute for Cosmic Ray Physics, Osaka City University, Osaka 558-8585, Japan 
2 KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
3 Department of Physics, Nagoya University, Nagoya 464-8601, Japan 
4 Department of Physics, Nara Women's University, Nara 630-8506, Japan 
5 Department of Physics, Tokyo Institute of Technology, Tokyo 152-8551, Japan 
6 Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 
7 Department of Physics, Purdue University, West Lafayette, IN 47907, USA 
8 Department of Applied Physics, Tokyo University of Agriculture and Technology, Tokyo 

184-8588, Japan 

9 Department of Physics, The Graduate School of Science and Technology, Kobe University, Kobe 

657-8501, Japan 

10 Department of Physics, Faculty of Liberal Arts, Tezukayama Gakuin University, Nara 631, 

Japan 

11 Department of Physics, Osaka City University, Osaka 558-8585, Japan 

Abstract 

We have studied inclusive muon events using all the data collected by the TOPAZ detector 
at y/s — 58GeV with an integrated luminosity of 273pb _1 . From 1328 inclusive muon events, 
we measured the ratio R q q of the cross section for qq production to the total hadronic cross 
section and forward-backward asymmetry A q FB for b and c quarks. The obtained results are 
R b l = 0.13 ± 0.02(stat) ± 0M(syst), R C c = 0.36 ± Om(stat) ± 0.05(syst), A b FB = -0.20 ± 
0A6(stat) ± O.Ol(syst) and A C FB = -0.17 ± 0A4(stat) ± 0.02(syst), in fair agreement with a 
prediction of the standard model. 
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1 Introduction 



The e + e _ — > qq cross section and charge asymmetry for heavy quarks (b and c quarks) are fun- 
damental quantities of electroweak interactions. Especially, in the TRISTAN energy region, the 
maximum forward-backward charge asymmetry for quark pair production, predicted by the stan- 
dard model, provides high sensitivity to quark couplings. In this paper, we report the final result 
of an analysis for the quark pair production cross-section ratio to the total hadronic cross-section 
(Rqq) and forward-backward asymmetry of b and c-quark (Afb), using the data of all the high- 
luminosity runs collected by the TOPAZ detector at the TRISTAN electron positron collider. The 
data were collected from 1990 to 1995 at y/s — 58 GeV, corresponding to a total integrated lumi- 
nosity of 273 pb _1 . Previously, by using parts of the TOPAZ data, we reported on measurements 
of the heavy quark production cross section and charge asymmetry: through inclusive muons(0), 
inclusive electrons (g [|) and D*±@). In the previous inclusive muon analysis, we derived the fa- 
quark parameters assuming the standard- model parameters for the c-quark with 41 pb _1 of data. 
The obtained result was A^ B = — 0.71±0.34(stat) ^Q Qg(syst)(^|). Through inclusive electrons, the 
results were Ap B = -0.49±0.20 (stat)±0.08(syst) and A^ B = -0.64±0.35(stat)±0.13(syst), using 
197 pb _1 of data(||). In this paper, we derive R q q and A-pb for both the b- and c-quark using the 
inclusive muon events with improved statistics. Data samples in previous analyses are included. 
In addition, we perform a correction on the pion punch-through rate using the measured data of 
pions from r pairs. 

The paper is structured as follows. A description of the TOPAZ detector and an overview of 
the data taking are given in Section 4 The hadronic event selection and the muon identification 
are described in Section |3|. Section ^ describes the analysis, including a Monte-Carlo simulation, 
flavor separation, fitting method, and systematic uncertainty of the measured parameters. Sec- 
tion |^ presents discussions of the fitting results and shows some figures compared to the other 
experimental results. Finally, conclusions are given in Section ^. 



2 Detector and data taking 
2.1 TOPAZ detector 

TOPAZ detector was a general-purpose Air detector located at the TSUKUBA experimental hall. 
A quadrant cross section of the TOPAZ detector is shown in Figure |l| The detector was upgraded 
in time for the high-luminosity runs which started in 1990. The upgrade was done by adding 
a vertex chamber, a ring calorimeter and forward-backward muon chambers; further, the inner 
drift chamber was replaced by a trigger chamber, and the luminosity monitor was replaced by a 
forward calorimeter. After the upgrade, tracking of charged particle was done by a vertex chamber 
(VTX) (|5]), a trigger chamber (TCH) and a time projection chamber (TPC)(||), which were placed 
inside of a IT magnetic field, produced by a super-conducting solenoid magnet (SCS) ((?]). Time- 
of-flight counters (TOF) (|) placed inside of SCS were used to provide information about the 
timc-of- flight and trigger. The energies of electrons and photons were detected with a barrel 
calorimeter (BCL)(|), a ring calorimeter (RCL), an endcap calorimeter (ECL)([l0|; |Tl| ) and a 
forward calorimeter (FCL) ([l2]) . Those calorimeters were installed outside of the tracking devices. 
The total angular coverage of these calorimeters was | cos#| < 0.998. The muon detection system 
(MDC)(|l^) consisted of the barrel part and the forward-backward part. They were placed at the 
outermost part of the TOPAZ detector. The coordinate system used was: z for the direction of the 
electron beam, 9 for the polar angle measured from the z axis, (f> for the azimuth angle measured 
from the horizontal direction pointing to the outward direction of the accelerator ring and r for 
the radial direction from the beam axis. In this analysis, we mainly used TPC, BCL, ECL and 
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Figure 1: Quadrant cross section of the TOPAZ detector. 



MDC. 

Time projection chamber (TPC) The Time Projection Chamber (TPC) was the central 
tracking device. Its function was to measure the 3-dimensional positions and dE/dX of charged 
particles. The inner and outer radii of the TPC were 36 cm and 109 cm, respectively, and its 
length was 260 cm. Its polar angle coverage was |cos0| < 0.83. They were made of 8 sectors 
divided into the forward and backward side, separated by a central membrane. The detection of 
drifted electrons was performed by 16 multiwire proportional chambers, placed on the end plane 
of each sector. The chamber was filled with an Ar/CIT = 90/10 mixture gas (P10) at 3.5 atm. 
The momentum resolution of TPC was o Pt /pt = \/(l-5pT(GeV/c)) 2 + (1.6) 2 %. 

Electromagnetic calorimeter (BCL and ECL) The barrel calorimeter (BCL) was made 
of 4,300 lead glass blocks with photo- multipliers. Each lead glass block had a 20 radiation 
length. Its polar-angle coverage was |cos#| < 0.84. The BCL energy resolution was ue/E = 
^(8.0) 2 /£(GeV) + (2.5) 2 %, measured by Bhabha events. 

The endcap calorimeter (ECL) was a pair of sampling gas calorimeters, and it consisted of a 
sandwich of conductive plastic tubes with cathode readouts and 2 or 3 mm thick lead plates. The 
total radiation length of the detector was 18. Its coverage was 0.85 < | cos#| < 0.98. The energy 
resolution of ECL was <je/E = 6.7% for Bhabha events. 
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Muon drift chamber (MDC) The muon detector (MDC) consisted of a barrel part (BMU) 
and a forward-backward part (EMU). In the present analysis, we only used BMU. BMU was 
located outside of the return yoke, which was made of 40 cm thick steel. BMU consisted of 3 
superlayers of muon chambers interleaved with 2 layers of 30 cm thick muon filters made of steel. 
Each superlayer had double-layers of drift chambers, both measuring the <fi direction of the track. 
In the outermost superlayer, we had an extra double-layer of drift chambers that measured the 
z direction of the track. The polar-angle coverage of BMU was | cos6*| < 0.66, and the minimum 
absorption length of BMU was 7.9. The muon drift chambers were made of extruded aluminum 
tube modules. One module contained eight cells, each 10 cm wide and 5 cm high, arranged in two 
layers of four cells each, staggered by a half cell to cover the dead regions near the cell walls. The 
averaged detection efficiency of each plane was 96%, measured using cosmic rays. 

Trigger system The main triggers were made of an energy trigger, a track trigger and a muon 
16|). The majority of the hadronic events were triggered by the energy trigger. The 
energy trigger required one of the following four conditions to be satisfied: 1) The total energy 
deposited in BCL be greater than 2 or 4 GeV, depending on the run conditions, or 2) the sum of 
the energy in the forward and backward ECL detectors be greater than 10 GeV, or 3) the BCL 
system had two energy clusters, both greater than 1 GeV, or 4) there were no energy clusters 
of energy greater than 3 GeV in neither the forward nor the backward ECL. The rate of energy 
triggering was approximately 1 Hz. 

The track trigger required that there be at least 2 tracks having an opening angle larger than 
45° in the r — <f> plane. Their vertex position was required to be within ± 20 cm of the interaction 
point along the beam direction by a software trigger. 

The muon trigger system was implemented for triggering e + e~ — > events in high- 

luminosity runs. Muon track signals were made of 8-OR signals of BMU drift tubes in coincidence 
with TOF signals. The muon trigger required two muon track signals. The opening angle of two 
tracks was required to be greater than 135°. 

2.2 Data taking 

TOPAZ started data taking in the spring of 1987 and ended in the summer of 1995. The total 
integrated luminosity was 340pb _1 . In 1990, the TRISTAN main ring was upgraded so as to 
increase the luminosity. Then, the collision energy was fixed at 58 GeV to maximize the luminosity. 
The data used in this analysis were taken at y/s = 58 GeV during the period given in Table U 
The data before 1990 and the data in January 1993 were not used because the beam energy was 
not 29 GeV. The data in the first quarter of 1995 were not used because of a detector problem. 



Table 1: Event sample. 



Period 




integrated luminosity (pb *) 


No. of Hadronic events 


Feb, 1990- 


-Dec, 1992 


114 


12,811 


Feb, 1993- 


-Dec, 1994 


144 


14,942 


Apr, 1995- 


-May, 1995 


15 


1,808 






273 


29,561 



3 Inclusive muon selection 



trigger (04; 15 
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3.1 Hadronic event selection 



To select inclusive muon events, we first selected hadronic events from DSTs using the TOPAZ 
standard selection criteria for hadronic events ( P~7[ ) , which were: 

(1) At least five "good" tracks were coming from the interaction point, where a "good" track 
was defined by (i) r < 5 cm and \z\ < 5 cm at the closest point of approach to the beam 
axis, (ii) p T > 0.15 GeV/c, and (iii) | cos#| < 0.83. 

(2) The total visible energy, E V i S , had to exceed the beam energy, E\> eam . 

(3) The momentum balance, | ^2p z \/ E V i Sl was less than 0.4. 

(4) The larger of the invariant jet masses in the two hemispheres, divided by the plane perpen- 
dicular to the thrust axis, Mj et , exceeded 2.5 GeV/c 2 . 

(5) The number of clusters having an energy greater than £?b eam /2 was less than 2. 

The background processes for the hadronic events were mainly r pairs and two-photon events. 
Especially, r pairs have a 17% branching ratio for decaying into muons, and have a large forward- 
backward asymmetry; hence, it could cause a significant background to the asymmetry measure- 
ment. To reduce these backgrounds, we required additional tight selection criteria, which were: 

• At least two charged tracks existed in each hemisphere. 

• A tighter cut on the invariant jet mass in the criteria(4): Mj e t > 3.5 GeV/c 2 . 

A total of 29,561 events passed the above requirements. Using Monte-Carlo simulations, we 
determined the efficiency of this hadronic event selection to be 65.08 ± 0.07%(systematic error 
only). The remaining background in the selected events were estimated to be 0.2% from r pairs 
and 0.1% from two-photon events, according to a Monte-Carlo study ( |lS| ; [l9| ). These are negligible, 
compared with the statistical errors and the other systematic errors in this analysis. Out of this 
hadronic event sample we excluded those events that had no high voltage on MDC. After this 
selection, 27,614 hadronic events remained. 

3.2 Muon identification 

Muon identification in hadronic events suffered from fake muons, such as punch-throughs and 
decay-in-flights of hadrons. In order to discriminate prompt muons from fake muons, following 
criteria were employed: 

(i) The track had to pass the "good" track cuts, described in the previous subsection. 

(ii) The momentum of the track had to be greater than 2.5 GeV/c and | cos#| < 0.6. 

After the track passed (i) and (ii), the track was extrapolated to MDC, assuming that it was 
a muon. We then applied the following. 

(iii) All three r — 4> superlayers of MDC had to have at least one hit in each superlayer within 20 
cm or S&track from the extrapolated tracks. Here, <jf rack is the quadratic sum of the track 
extrapolation errors, the multiple Coulomb scattering errors, detector space resolution, and 
the error due to the detector alignment. In addition, we required that at least one superlayer 
had to have adjacent firing cells in the other layer of the same superlayer. 

(iv) MDC hits were not allowed to be shared by the other tracks. 
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Figure 2: Ad (a) and Ad/a tra ck (b) for muons. The dots with error bars are cosmic-ray data. 
Histograms are Monte Carlo prediction. 



The capability of the Monte-Carlo to simulate muons has been tested with cosmic-ray muons. 
Figure || shows the distance between the extrapolated track and the associated hit, Ad, and 
Ad/ a 'track- As can be seen in these figures, the Monte-Carlo simulation and the data are in 
good agreement. The muon identification efficiencies as a function of momentum p and also as a 
function of cos 6 are shown in Figure [| The efficiency shows a plateau over 2.5 GeV/c in p and 
| cos#| < 0.6. The measured efficiency for muon identification is 93% at the plateau, which agrees 
with the Monte-Carlo prediction within 1%. The inefficiency was due to dead wires in MDC (3%), 
dead regions in the sector boundaries (3%) and cuts on the distance between the MDC hits and 
the track (1%). 

The number of inclusive muon events was 1,328, after selecting the hadronic events by muon 
identification. If an event had two or more muon candidates, we chose the one with the highest 
momentum as the muon for tagging c and b quarks. 



4 Analysis 

4.1 Monte Carlo simulation 



We used JETSET7.3(|20|) for e+e" -> qq event generations. The used parameters in JETSET7.3 
were tuned-up for hadronic events using the event shape data (plf). In the hadronization, we 
used the LUND symmetric function for the fragmentation function for light quarks (u,d,s) with 
a = 0.413 and b — 0.9. For heavy quarks (c,b), a function by Peterson et a/.(|22]) was used with 
e c = 0.05 and e& = 0.01. For the standard model parameters, we used sin 2 0\y = 0.2315, M^o = 



91.187 GeV/c 2 , and T z0 = 2.490 GeV(23) 



For a detector simulation, we used the TOPAZ detector simulator, which simulated the be- 
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Figure 3: Muon identification efficiency measured through cosmic-ray muons (dot with error bars) 
and Monte-Carlo prediction (dotted line) as a function of the momentum (a) and polar angle (b). 



haviors of the particles in the TOPAZ detector: such as the energy loss, multiple scattering, 
decay-in-flights and the detector signals. For the simulations of particle's interactions with the de- 
tector material, EGS4(p4|) was used for electromagnetic processes and GHEISHA-7(p5|) for nuclear 
interactions. We used 314,463 hadronic events for the studies described in this paper. 

To test the validity of the hadronic event selection and the Monte-Carlo simulation, we exam- 
ined the general features of hadronic events. Figure |] shows the distributions of the momentum (a) 
and the polar angle (b) of the track, and the polar angle of the thrust axis (c) . The data and the 
Monte-Carlo results are in good agreement, except for a dip observed at 9 = 90° in the polar-angle 
distribution. This dip is due to the effect of the central membrane of TPC. The track reconstruc- 
tion inefficiency at the membrane is 0.4% for hadronic events and 0.8% for inclusive muon events. 
The effect of this inefficiency is negligible on the cross section and asymmetry results. 



4.2 Background estimation 

Possible background sources to prompt muons are hadron punch-throughs, muons from decay-in- 
flights of light hadrons (mainly tt^ and K ± ) and accidental hits in MDC due to the beam back- 
grounds or cosmic rays. To estimate the rate of accidental hits, we applied the muon identification 
criteria to the electron (positron) tracks in Bhabha events, obtained in the same experimental pe- 
riod. No track was identified as a muon. From this result, the accidental hit rate was considered 
to be negligible. 

The other background sources were evaluated using Monte-Carlo simulations. Since the pion 
and kaon lifetimes and decay modes are well established, the decay-in-flights of light hadrons 
can be precisely calculated in the simulation. On the other hand, the rate of hadron punch- 
throughs depends on nuclear interactions with detector materials. Therefore, it must be verified 
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Figure 4: Hadronic events; momentum (a) and polar angle (b) of good tracks, and angular distri- 
bution of thrust axis (c). The dots with error bars are the data; the histograms are Monte-Carlo 
results. 
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Figure 5: Punch- through rate as a function of momentum derived through the pions from r-pairs. 
The dots with error bars are the data. The dotted and dashed lines are Monte Carlo simulations 
before and after tuning the parameters in GHEISHA. 



experimentally. For this purpose, we studied the punch-through of pions by using 1,202 charged 
pions from r — ► mnrv and r — > pv — > tt^tt ^ in the e + e~ — * t + t~ events. Figure ^ shows the 
experimental ratio together with the Monte-Carlo expectation as a function of the momentum. 
As can be clearly seen in this figure, it was found that the Monte-Carlo prediction significantly 
underestimated the rate, especially in the high- momentum region with the default parameters. To 
correct this difference, we tuned up the most influential parameter, which is the total cross section 
of pions (kaons) interacting with materials in the GHEISHA routine. We introduced a correction 
factor to scale the total cross section as a function of the momentum, as listed in Table |2[ In the 
high-momentum region of p > 10 GeV/c, the correction factor was determined through the pion 
punch-through rate in the r pair events. For p < 3 GeV/c, the correction factor was estimated 
through energy deposits of hadrons in BCL, and it was fitted to a function. These two correction 
factors were linearly connected in the intermediate momentum regions, i.e. 3 < p < 10 GeV/c. 

The results of the tuned-up Monte-Carlo simulation are also shown in Figure |. The discrep- 
ancy between the Monte-Carlo results and the data was reduced significantly with the tuncd-up 
Monte-Carlo procedure. 



Table 2: Correction factor for the cross section for each momentum region. 



Momentum region (GeV/c) 


correction factor f(p) 


p < 3 


f{p) = -0.013^ + 0.139p + 0.826 


3 < p < 10 


f(p) = -0.054p+ 1.295 


10 <p 


f( P ) = 0.726 
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Figure 6: Distributions of momentum (a) and polar angle (b) of inclusive muons. The dots with 
error bars are data and the histograms are Monte-Carlo results. The open regions show the 
contributions from c and b quarks. The shaded regions are the background from light quarks. 



In the inclusive muon sample, the fractions of decay-in-flight and hadron punch-through were 
estimated to be 25% and 23%, respectively, by using the tuned-up Monte-Carlo program. Figure 
^ shows the momentum spectrum (a) and polar-angle distribution (b) . The background from light 
quarks was obtained by a Monte-Carlo simulation using the total number of hadronic events for 
normalization. 



4.3 Flavor separation 

Inclusive muon events can be categorized into four groups by the sources: (A) prompt muons 
from primary b-decays (b — > fi), (B) cascade muons from b — > c(c) decays (b — > c — > fi), (C) 
prompt muons from primary c-decays (c — > fi), and (D) backgrounds from decay-in- flights and 
punch-throughs (others). In this analysis, we included the b — > t — > \i mode in category (A). 
Table || gives the results of Monte-Carlo studies, giving the fraction of each components in the 
inclusive muon events. 

We used the thrust axis to determine the primary quark direction. The accuracy of the quark 
direction determined by this method was estimated to be 5°, based on a study by a Monte-Carlo 
simulation. We define the angle of the quark direction by — Q cos8 trus t, where Q is the charge 
of the tagged muon and 9 is the angle of the thrust axis with respect to the beam axis. To 
separate the quark flavors, we used the muon transverse momentum with respect to the jet axis, 
Pt*- Jets were reconstructed using the JADE jet clustering algorithm (g6|) with a scaled invariant 
mass cut, Y cut (= Mij/E vis ) = 0.04. Muons from b-quarks have a larger p^f compared with those 
from c-quarks, due to the heavy b-quark mass. Using the p 3 ^ cut at 0.8 GeV/c, we classified 
the inclusive muon events into a b-enriched sample (p^ > 0.8 GeV/c) and a c-enriched sample 
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Figure 7: p^T* distribution of inclusive muons (a) and purities of b and c-quarks as a function of 
Py* cut (b). In (a), dots with errors show the data. Histograms are Monte Carlo: b — ► p(open), 
b — ► c — > //(hatched), c — > //(double-hatched) and others (triple- hatched), respectively. In (b), 
solid curve shows the b-quark purity and dotted curve shows the c-quark purity. 



(p^ < 0.8 GeV/c). We chose p^r* = 0.8 GeV/c, because the systematic uncertainty was the 
smallest at 0.8 GeV/c. The purity of the c- and b-quarks as a function of the pr cut was studied 
by a Monte-Carlo simulation. Figure ^ shows the b, c-quark purities as a function of the px cut. 
For b-quarks, the purity was calculated by summing up all of the b-quarks with px greater than 
a given pt cut, divided by the total number of inclusive muon tracks in the same pt region. For 
c-quarks, summing was done for px less than a given px cut. 



Table 3: Percentage of muons from each source in the inclusive muon sample. 



Source 




Fraction(%) 




p> T et < O.f 


5 GeV/c y T et > 0.8 GeV/c 


total 


b — > jU 


9.7 


41.9 


20.1 


b — > c — > /i 


5.8 


5.2 


5.6 


C — > jU 


31.6 


13.2 


25.7 


others 


52.8 


39.6 


48.5 



4.4 Simultaneous fit for b and c-quark 

In a previous paper (|l|), we derived only the b-quark parameters by assuming the standard- model 
parameters for the c-quark because of low statistics. This time, we performed a four-parameter fit 
of (R b i, R C c, A b FB , A C FB ) top and —Q cos 6 thrust ipr — 0.8 GeV/c and < 0.8 GeV/c) distributions. 
Definitions of R q q, and Apb are given by 

Rqq = a(e + e~ — ► qq)/a(e + e~ — > hadrons) (1) 
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and 

A FB = ^-^ (2) 

where a represents the lowest order (SU(3) x SU(2)l x U(l)) cross-section and (Jf{&b) is the 
cross-section for quark to travel forward (backward) with respect to the e~ direction. The \ 2 of 
the fit was defined by 

2 _ 2 , 2 , 2 

Xtotai — X C ose(p^ et >0.8Gey/c) ^cos6((p^ !t <0.8Gey/c) 



»=i < 3=1 ^; k=i a K 



(3) 



where iV,- and a N ^ are the number of inclusive muon events and the statistical error for each bin, 
respectively. is the number of expected events for each bin, which is defined by 

K = NZ P d {R b - b ^Br^(b^ f ,)F(A b FB ) l -Ct t2 

+ R bl ■ 2Br(c - M )(l + (_i^^ B ) 4 ■ C*^~" 

+ i? ce • 2Br(c -> fj,)F(-A c FB )i • C,r M } 

+ pothers ( 4 ) 



for -QcosOthrust bins, and 



+ i? h5 ■ 2Br(c -» M )(l + a) ■ W^ c ^ 
+ i? ce ■ 2Sr(c -» M ) • W^} 

+ ]\[Others ^ 

for p bins, where iVf^? is the total number of hadronic events. Br sum (b — > /i) is the combined 
branching ratio of i?r(6 — > /i) and Br(b — > r — ► /i). £?r(fo — > /i) and i?r(c — * /j) are the branching 
ratios of 6 — > /i and c — * /x, which are equivalent to the branching ratios for muons semileptonically 
decayed from b and c hadrons, respectively. Br(b — > ccs) is denoted as a in the above equations. 
The used branching ratios are listed in Table |[ F(A FB )i is the polar-angle distribution function, 
integrated in each cos# bin (i-th bin), which is written as 

F(4 B )i= /|(l + cos 2 0+|^ B cos0)dcos0. (6) 
Ji ° -J 

(jmode j g correction factor for the each cos 9 bin from the decay mode, given by 

C1 i x\iq 

°* — fl i X\had „had £ i ' V ' J 



where (1 + <5)* 9 is the combined correction factor of the QED (initial state photon radiation) and 
QCD (parton shower) radiative corrections for the i-th bin of e + e~ — > gg. Similarly, (1 + 8)^f al 
is a radiative correction factor for the total hadronic cross section. These correction factors were 
estimated through a Monte-Carlo simulation based on the LUND event generator, JETSET7.3. 
r\\ q is the acceptance correction for the i-th bin of each mode, rj had is that for hadronic events, 
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Table 4: Branching ratios used in the fits. 



Mode 


Branching ratio(%) 


Br sum (b -> u) 


11.3 ±Q.5 


Br(b -> m)@ 


10.8 ±0.5 


Br (6 r-> M)@ 


0.45 




9.0 ±0.7 


-Br (6 — > ccs)(|30|) 


21.9 ±3.7 



and is the muon identification efficiency for the i-th bin (cos# bin). \y™ ode is the correction 
factor for each p bin, which is given by 

-ijrmode _ (1 ± ^Yk ^k,mode u , Q \ 



Table | gives the correction factors for each bin. N° thers is the number of background events 
for the i-th bin, which is also listed in Table |5|. The numbers were derived through Monte Carlo 
simulations. 

As a result of the four-parameter fit (Figure || and ||), we obtained R b i = 0.13 ± 0.02, R cS — 
0.36 ±0.05, A b FB = -0.20 ±0.16 and A C FB = -0.17 ±0.14, with x 2 /D.O.F = 14.93/18. The errors 
are statistical only. The correlation coefficients obtained from the fit is given in Table §. The 1 
contours of the fit for R cc v.s. Rbb and A c v.s. At, are shown in Figure ^|. 

4.5 Systematic errors for R q q and A q FB 

Possible sources of systematic errors are listed in Table 0. The largest systematic error comes 
from the uncertainty in the probabilities of mis-identified hadrons as muons. This is caused by the 
light quark background in the heavy quark samples. In order to estimate the effect of hadron mis- 
identifications, we performed the same analysis by changing the correction factor of the effective 
pion- nucleus cross section by ±lcr of the measured error of fake rate. The effect of the cuts on 
muon identification was studied by varying the cut values on momentum p, Ad, and lS.d/<r t rack 
by ±0.5 GeV/c, ±5 cm, and ±0.5er trac /j, respectively, and examined the change in the obtained 
cross sections and asymmetries. We studied the effect of the p 3 ^ cut by shifting the cut value 
by ±0.1 GeV/c. The effect of the uncertainly in the MDC detector acceptance was checked by 
varying the | cos tra ck \ cut from 0.6 to 0.58. Based on these studies, we selected the cut values 
that minimize the systematic errors. The error due to the semileptonic branching ratios for b and 
c hadrons were obtained by changing the branching ratios by ±1<t of the quoted number. Changes 
of fragmentation parameters cause the changes in the muon identification efficiency and flavor 
separation. The error due to this effect was estimated by changing the fragmentation parameters 
by ±lcr of the quoted number. 
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Figure 8: Results of the four-parameter fit: cosOthrust for the c-enriched sample(p^T* < 0.8 GeV/c) 
(a) and the b-enriched sample(p^f > 0.8 GeV/c) (b), and momentum distributions (c). The dots 
with error bars are the data. The histograms are the results from the four-parameter fit: b — > 
/x(open), b — > c — > ,u(hatched), c — > ^(double-hatched) and others (triple-hatched), respectively. 
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Figure 9: Results of the four-parameter fit for R q q (a) and A q FB (b). The cross bar indicates the 
central value and the ellipse shows la contour. The errors are only statistical. The points show 
the standard-model prediction. 
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Tabic 5: Number of events and correction factors in each bin. 



(a) 


Correction factor for each cos 


6 bin: 


p^ et < 0.8 


(GcV/c) 


cos 


6 bin 


Nt 


Ajothers 


c b^ 




CT" 


-0.8 


- -0.6 


20.± 4.5 


14.0±1.1 


0.049 


0.034 


0.038 


-0.6 


- -0.4 


133.±11.5 


65.0±2.4 


0.178 


0.180 


0.271 


-0.4 


- -0.2 


142.±11.9 


76.2±2.6 


0.222 


0.187 


0.319 


-0.2 


- 0.0 


152.±12.3 


76.9±2.6 


0.183 


0.137 


0.288 


0.0 


- 0.2 


154.±12.4 


76.6±2.6 


0.234 


0.149 


0.273 


0.2 


- 0.4 


132.±11.5 


76.5±2.6 


0.201 


0.176 


0.313 


0.4 


- 0.6 


144.±12.0 


60.5±2.3 


0.209 


0.128 


0.253 


0.6 


- 0.8 


24.± 4.9 


9.8±0.9 


0.064 


0.029 


0.044 



(b) Correction factor for each cos 6 bin: p^ et > 0.8 (GeV/c) 



cos 


6 bin 


N? 


j^others 








-0.8 


- -0.6 


22.±4.7 


4.8±0.7 


0.106 


0.042 


0.014 


-0.6 


- -0.4 


62.±7.9 


22.0±1.4 


0.326 


0.062 


0.046 


-0.4 


- -0.2 


80.±8.9 


26.3±1.5 


0.449 


0.057 


0.070 


-0.2 


- 0.0 


64.±8.0 


24.8±1.5 


0.476 


0.063 


0.062 


0.0 


- 0.2 


64.±8.0 


26.3±1.5 


0.443 


0.073 


0.050 


0.2 


- 0.4 


66.±8.1 


25.2±1.5 


0.411 


0.080 


0.060 


0.4 


- 0.6 


54.±7.3 


17.4±1.2 


0.334 


0.037 


0.040 


0.6 


- 0.8 


15.±3.9 


5.7±0.7 


0.125 


0.014 


0.008 



(c) Correction factor for each momentum bin 



Momentum bin (GeV/c) 


K 


j\jothers 








2.5 -3.0 


181.± 13.5 


95.0 ±2.9 


3.048 


2.902 


2.998 


3.0 -4.0 


281. ± 16.7 


116.9±3.2 


4.921 


3.764 


4.903 


4.0 -6.0 


325. ± 18.0 


124.6±3.3 


4.572 


1.573 


2.627 


6.0 -10.0 


309. ± 17.6 


130.5±3.4 


2.600 


0.402 


0.936 


10.0-16.0 


165. ± 12.8 


98.2 ±2.9 


0.865 


0.045 


0.145 


16.0-30.0 


68. ± 8.2 


43.2 ±1.9 


0.091 


0.007 


0.017 



Table 6: The correlation coefficients for the parameters in the four-parameter fit. 





Rbb 


Rcc 


^FB 


A c 

FB 


Rbb 


1 


-0.77 


0.01 


-0.17 


Rcc 




1 


-0.37 


0.16 


A b 

FB 






1 


0.44 


^FB 








1 
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Table 7: Summary of systematic errors 



Source 




AR cZ /R c5 


&A» FB /A FB 


AA C FB /A FB 


hadron mis-ID 


4.9% 


12.1% 


3.8% 


13.3% 


muon ID 


<0.1% 


0.6% 


1.5% 


0.9% 


Py* cut 


0.9% 


0.6% 


1.5% 


0.2% 


MDC acceptance 


<0.1% 


<0.1% 


0.5% 


0.9% 


branching ratio 


2.8% 


5.0% 


2.1% 


1.2% 


fragmentation parameter 


1.7% 


3.8% 


<0.1% 


<0.1% 


total 


6.1% 


13.7% 


4.9% 


13.8% 
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5 Discussion 



The results of the four-parameter fits with systematic errors are R bb = 0.13±0.02(siai)±0. 01 (syst), 
R cS = 0.36 ± 0M(stat) ± 0.05(syst), A b FB = -0.20 ± 0.1Q(stat) ± O.Ol(syst) and A C FB = -0.17 ± 
0.14:(stat) ± 0.02(syst). Four-parameter fits of the c- and b-quark have correlations in the param- 
eters of the c-quark and b-quark. For example, an increase of R c£ would cause R b i to decrease. 
Also, an increase of A C FB would cause A b FB to increase. To estimate this effect, we applied one 
quark parameter fit, by fixing the other quark parameters to the standard-model predictions. The 
used values of the standard model parameters are R bb = 0.13, A b FB = —0.43 for the c-quark fit and 
R cS = 0.30, A C FB = -0.48 for the b-quark fit. For A b FB , correction for B° d s B° d s mixing is included 
using the mixing parameter at high energy, \b = 12%. The results of the fits are R bb = 0.15±0.01, 
A b FB = -0.29 ± 0.13(b-quark fit) and R cS = 0.36 ± 0.03, A C FB = -0.26 ± 0.13 (c-quark fit). The 
results are summarized in Table 0. They are consistent with the four-parameter fit. 



Table 8: Summary of R q q and A q F 



fitting method 


flavor 


Rqq 


^FB 


X VD.O.F 


four-parameter fit 


c-quark 
b-quark 


0.36 ±0.05 
0.13 ±0.02 


-0.17±0.14 
-0.20 ±0.16 


14.93/18 


1 quark parameter fit 


c-quark 
b-quark 


0.36 ±0.03 
0.15 ±0.01 


-0.26 ±0.13 
-0.29 ±0.13 


19.29/20 
19.29/20 


Standard Model predictions 


c-quark 
b-quark 


0.30 
0.13 


-0.48 

-0.43 (mixing X b = 12%) 





The errors are statistical only. 



The differential cross sections for b-quark production (d<j bb /d cos 9) and c-quark production 
(da eg / 'd cos 9) were obtained from the cos Othrust distribution. In order to derive da bb /dcos9, the 
contributions from light quarks and cc were subtracted, and bb was subtracted for da c£ /dcos9. 
The results are given in Table |^ and Figure [w| The measured asymmetry for the c- and b-quarks 
are smaller than the standard-model prediction by 2.2<r and 1.5<r, respectively. 

Figure [ll] shows the measured R cS and R bb together with the other experimental dataj27j; 
|3i|; |32|). The measured asymmetries of b and c-quark are shown in Figure O with the other 



experimental data(§ [§ ||; ||; |J; HI 




1 □ B 




9; 50: 51) 



TGTlj 



The asymmetries measured at PEP(|36|; ^jj |38| ^0] ^2(), PETRA(L_ 

|5(j| : pl| ), and LEP(|27|) are combined for each experiment. We combined the present results of A q FB 
with our previous measurement by inclusive electron ( for c and b quark ) and D*^ (for c quark). 
The combined values are A b FB — —0.28 ± 0.15, A C FB — —0.35 ± 0.09. The combined results are 
also shown in Figure n3. 



6 Conclusion 

We have studied inclusive muon events and extracted the cross-sections and charge asymmetries 
of b and c-quark using all of the data collected with the TOPAZ detector at = 58 GeV. A 
total of 1,328 inclusive muon events were selected from 29,561 hadronic events with an integrated 
luminosity of 273 pb _1 . To improve the accuracy of the simulation for pion punch-throughs, we 
tuned the parameters in the Monte-Carlo simulation, using 1,202 pions from — > 7r ± 7r ± 7r =F ^, and 
— > p ± i i 1 p ± — > tt^tt in e + e~ — » t + t~ reaction. The measured ratio, Rqq, of the cross section 
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cose 



Figure 10: Differential cross section of e + e~ — ► bb and e + e~ — > cc (upper figures). The solid lines 
are the best fits and the dashed lines are the standard model predictions. Lower figures show the 
difference between the observed cross section and the standard model prediction. The open circles 
are the result of the four-parameter fit and the cross marks are the result of the 1-quark parameter 
fit. 
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Figure 11: R cS (a) and R b i (b) as a function of the center-of-mass energy. 
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Figure 12: Forward-backward charge asymmetry for c-quark (a) and b-quark (b) as a function of 
the center-of-mass energy. 
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Tabic 9: Differential cross sections for e + e — > cc and bb processes. 



cos 


bin 


1 da cE 


i d °bh 


<Jh n d COS 9 


Ohnr] d COS 


-0.8 ~ 


-0.6 


0.153±0.113 


0.088±0.036 


-0.6 ~ 


-0.4 


0.238±0.049 


0.073±0.020 


-0.4 ~ 


-0.2 


0.101±0.038 


0.074±0.016 


-0.2 - 


0.0 


0.200±0.047 


0.046±0.014 


0.0 ~ 


0.2 


0.192±0.044 


0.047±0.015 


0.2 ~ 


0.4 


0.119±0.036 


0.038±0.016 


0.4 ~ 


0.6 


0.148±0.044 


0.062±0.018 


0.6 ~ 


0.8 


-0.057±0.122 


0.037±0.025 



for qq production to the total hadronic cross section and the forward-backward asymmetry, A q FB , 
of b and c quark arc R bh = 0.13 ± 0.02(stat) ± O.Ol(syst), R c5 = 0.36 ± 0.05(stat) ± 0.05(syst), 
A b FB = -0.20 ± O.W(stat) ± O.Ol(syst) and A C FB = -0.17 ± O.U(stat) ± 0.02{syst), respectively. 
The standard-model prediction for those parameters with 12% B — B mixing (xb) are R^i = 0.13, 
R C c — 0.30, A b FB = — 0A3,A C FB = —0.48, respectively. The measured R b i and R c5 agree with the 
standard-model predictions, though the measured A b FB and A FB are smaller by 2.2er and 1.5a, 
respectively. 
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